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Evaluation of Turbulence Induced Noise in Coherent Anti-Stokes 


Raman Scattering* 


ABSTRACT 

The effect of turbulence in & transonic wind tunnel on 
coherent anti-Stokes Raman scattering is considered. The driving 
pump and Stokes waves are taken to be coaxially propagating 
Gaussian beam waves which are focused on the Raman active medium 
through the turbulent boundary layer of the flow tube. The random 
index of refraction variations in the layer are modelled as phase 
perturbations of the driving waves which cause a reduction of the 
mean on-axis field and an increase in the mean diameter of the 
beams. Effective Gaussian beam parameters are developed and the 
radiated anti-Stokes power calculated as a function of the phase 
screen parameters. It is found that a significant reduction in 
signal strength occurs for realistic estimates of the phase screen 
parameter appropriate to a confined transonic flow. A method for 
estimating the signal degradation which could be applied to other 
experimental situations is presented. 

*This work has been submitted for publication to the Journal of 
the Optical Society of America under the title ’’CARS Through a 
Turbulent Boundary." 


I, Introduction 

•) I ■■■ lai BWaMnnHMMMMNM 

This project was conc^(rned with the effects of turbulence on 
the coherent anti-Stokes I’laman scattering (CARS) p'S'ocesa as it 
might be used to monitor conditions in a transonic wind tunnel. 
CARS has already proven to be a valuable tool for remote, 
non~intrusive, sensing of molecular species concentration and 
temperature in flames,^ combustion chambers^ and high speed 
flows,^»5 The major reasons for the great interest in CARS as a 
system probe stem from the facts that it does not perturb the 
system; it is capable of providing both spatial and temporal 
resolution; the signal is coherent and unidirectional thus 
providing excellent discrimination against background radiation; 
and both temperature and number density information can be 
extracted from it simultaneously. 

It has be;?n recently suggested that because of these 
properties CARS might be used to monitor conditions in transonic 
wind tuvmels near test air foils. Since the gas flow under such 
circumstances is often turbulent, the effect of the random index 
of refraction variations due to turbulence on the CARS process is 
a matter of some importance. Since the flow in wind tunnels 
is most turbulent near the walls of the flow tube the situation 
can be modelled as a thin turbulent layer r^iur, bounding a quiescent 
laminar flow region. 

CARS may be described in semi-classical macroscopic terms as 
a degenerate four wave mixing process mediated by a nonlinear 
susceptibility.^ Two waves, a pump wave at frequency and a 
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Stokes wave at ftequency ^2 shifted from the pump by the frequency 
of a Raman transition - 0)2 •* v) drive the system and. induce a 
polarization at W3 » Wi + v, the anti-Stokes frequency. This 
polarization field results in an emitted radiation field with each 
volume element in the interaction region producing a field at 
frequency W3 with a particular phase. For certain propagation 
directions, namely those for which the wave vectors satisfy the 
equation 2ki - ka - ka " 0, these radiated fields add j.n phase to 
produce a coherent anti-Stokes signal. In gaseous systems, with 
inconsequential dispersion, the three waves may be colinear and 
still satisfy the phase matching condition although other 
arrangements such as BOXCARS^ are preferred in circumstances 
where spatial separation of the signal and driving waves is 
desired . 

The theory of CARS generation in a homogeneous medium when 
the driving waves are coaxial Gaussian beam waves has been 
developed by Bjorklund.® His work takes into account explicitly 
the effects of dispersion and focusing of the driving beams. More 
recently Guha and Falk^ have considered the problem from the 
point of view of optimal focusing and phase matching. The power 
of the anti-Stokes signal depends on the power of the Stokes beam 
and the square of the pump beam power. It increases 
asymptotically as the beams are more tightly focused and thus some 
loss of signal power is expected if for any reason the quality of 
the driving beams is i^graded. 
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A medium with a time varying index of refraction causes a 
propagating laser beam to be deflected in a random manner and its 
instantaneous diameter to be broadened, If the refractive index 
variations are severe enough the irradiance pattern can even be 
broken into several independent patches, A turbulent layer could 
thus be expected to have a profound effect on CARS, 

The problem of optical propagation in random media such as 
the turbulent atmospheric boundary layer has been thoroughly 
researched over the past two decades. One of the more fruitful 
formulations of the problem is the v?xtended Huygens-Fresnel 
theory^*^>^^ which allows one to describe the statistical moments 
of the optical field in terras of the spatial spectrum of the index 
of refraction fluctuations. 

The approach taken in the present work is to apply the 
Huygens-Fresnel formulation to CARS generation by coaxial Gaussian 
beams. First a brief review of the CARS process is presented and 
an explicit expression for the radiated anti-Stokes power 
developed. The effect of a turbulent layer on a propagating 
optical beam is then described and an expression for the 
anti-Stokes power derived when the pump and Stokes beams are 
perturbed by such a layer. Finally the loss in the CARS signal 
due to realistic estimates of the turbulence parameters in a 
transonic wind tunnel is determined. 
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II. CARS with Coaxial Gaussian Beam Waves 
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The development of the macroscopic theory of CARS presented 
here follows that published by Bjorklund.® The driving pump and 
Stokes waves at frequencies Wj, ( 1)2 are assumed to be lowest order 
Gaussian modes propagating coaxially along the z axis with 
identical confocal parameters and beam waist locations. The beams 
are assumed Co be linearly polarized in the same direction so that 
they may be treated as scalars and the total field written as a 
simple sum 

E(r,t) » Re[Ei(r) expC-i^it) + E 2 (r) exp(-iW 2 t)] . 

In the cylindrical coordinate system, transverse coordinate 
p, axial coordinate z, the pump and Stokes field amplitudes take 

1. IKC 4.V/LIU 

Ej(r) H Ej(p,z) 

= (k.I./nb)^'^^ (1 + 2iz/b)"^ exp[ik.z - k.p^/(b + 2iz>] 

J J J J 

( 1 ) 

2 

with the confocal parameter b = k;Wo- , Wo. being the 

J J 

1/2 

beam radius at the waist (z = 0). The factor (kjlj/b) 
normalizes the field so that the energy transmitted across any 
transverse plane 

00 4. 

2H / dPP E,(p,z) E. (P,z) = I., 

- J J J 


a constant, independent of how tightly the beams are focused. 
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The effect of the pump and Stokes fields in the nonlinear 
medium is to produce a polarization wave at a frequency 
0)3 » 20)^ - 0)2, i.e., 

^ -i^st 

P(r,t) » Re[P(r)e ] with 

P(r) “ 3/4 N 2 wi, - 0)3) Ei^(r) E2*(r) ( 2 ) 

where N is the number density of active atoms and X is tbe 
nonlinear susceptibility per atom. The radiation field EaCr) 
generated from this polarization wave may be obtained by 
decomposing P(r) into plane wave components with amplitude 

P(K) » / d^r P(r) exp(-iK*r). ( 3 ) 

Each polarization plane wave component generates a radiating plane 

wavei2 

E3(K;r ) = ( 2 iik ^/ka) {l - exp[-i(z^ + L) <K - ks + K2/2k3)]} 

0 0 o z 

X (K ' k3 + K^/2k3>“^ P(ic) expCi K*?^) ( 4 ) 

z o 

with k being the Vacuum wave vector and K = (k,K ). 
o • z 

The total generated radiation field is the sum of these 
components 

E3(r ) = / d^K E(K;r ). ( 5 ) 

0 o 

Using Eqs.( 2 )i ( 3 ) and ( 4 ) in Eq.( 5 ) and noting that the factor 
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{l - exp[“i(z^ + L)(k - U3 + t:V 2 k 3 )KK - ks + KV2k3)"^ 
OS? z 

can be written; 

"o 

/ exp{i[K - k 3 + K^/ 2 k 3 ](z’ - z )dz' 

-L ^ ° 

the result is 

E 3 (r^) * i 3/4 NX(2 it) ^(k^^/ks)/ d^K / ^ dz' / d^r exp [ik* (r^-r) ] 

•"L 

X Ei^(?) E 2 *(r) exp{i[K - ks + K^/ 2 k 3 ](z' - z )} (6) 

z o 

In this form the integral over may be performed immediately 
to give a factor 6(z' ~ z) . The z' integration then gives 

E3(?^) H E3(P^,z^) 

« i(3NXk^^)/(16 TT^g) / dz / d^p Ei^(r)E 2 (?) 

•*L 

X / d^K exp{-ik 3 '(z ~ z ) + i(z - z )K^/2k3 + xk*(P - P) } 

o o o 

(7) 

In gaseous systems such as flames or wind tunnels dispersion 
is so slight that coaxial beams may have essentially perfect 
wavevector matching so that ks “ 2kj - k2« When Ei(r) and E 2 (r) 
as given by Eq,(l) are substituted into Eq.(7) then the integrals 
over K and p may be evaluated to give 






0R!Q!s^nL m m tg 
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E3(r^) - i(3NXk^^)(4k3)"^ expCiksz^) 

z 

0 

X j dz (1 + 2iz/b)"*^ [2ki + ka ~ 2 ik 3 z/b]”^ [h(z ,z)]”^ 
-L ° 

X exp[-p^^/bH(z^,z) ] (8) 


where 


H(z iz) 
o 


1 -f 4z^/b^‘ 

2k 1 + ka “ 2ik3z7b 


2i 


(z - z ) 
o 

bk3 


(9) 


Equation (8) gives the radiated anti-Stokes field for any r^. The 
quantity observed in experiments is of course the irradiance and 
the quantity of interest is the total power crossing the exit 
plane, z » L_, of the system, i.e, 

t: 


I3(V - 


2TT / 


o o 





ao) 


Equation (10) along with (8) and (9) then provides the means 
to determine the strength of th? anti-Stokes signal. The effect 
of tighter focusing, length of the interaction region, etc., may 
be examined and the configuration of the experiment optimized. 

For purposes of this report one important fact emerges: almost 

all the anti-Stokes signal is generated over a region within a 
distance of 5 times the confocal parameter of the beam waists, 

-5b < z < 5b. This reduces the domain over which numerical 
integrations must be performed. 
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The optical field at an observation point in a homogeneous 
medium due to a finite source can be expressed in terms of the 
Kirchoff Integral which is the moiShematical statement of the 
Huygens-Fresnel principle. The field amplitude at 
r = (P,z)> i.e,, E(p,z), due to a source distribution in the plane 
e 5 ( 0 , 21 ), i.e., E(o, 2 i) with zi < z is given by*^ 

E(p,z) “ (k/2iri) COS'}'' / d^o jr - sj ^ exp(ikjr “ s ) E^(o,zi) (11) 

where ^ is the angle between r and the z axis. When the lateral 
dimensions of the source and the region of interest in the 
observation plane are much less than the distance from the source 
to the observation point the Fresnel approximations hold and 

E(p,z) “ [k/2Tri(z“zi) ]exp [ik(z-z 1 ) ] / d ^0 exp [ik| P-oj ^/2(z-zi) ] 

X E^(a,zi) . (12) 

The source distribution corresponding to a Gaussian beam wave 
in the plane z ® ~L which converging to a focus at z 0 is, from 
Eq.(l), 

E^(0,-L) = (kl^/TTb)^'^^ (l-2iL/b)"^ exp [-ikL-ka^/(b-2iL) ] . (13) 

If this expression for the source field is substituted into 
Eq,(12) and the integration performed Eq.(l) is recovered 
demonstrating the fact that the formulation is consistent. 

The effect of a thin layer of turbulent gas on a laser beam 
propagating through it as illustrated in Figure 1 is to introduce 
random phase changes across the diameter of the beam which then 


8 


ORIQINAL 
OF POOR QUAUTV 


affect both the amplitude and phase of thu field st locations 
remote from the turbulent layer. This situairion may be modelled 
by introducing a random phase perturbation in the source field in 
Eq,(12) with statistical properties appropriate to the index of 
refraction fluctuations in tne turbulent layer^ the source 

field to be used in Eq,(12) is written as 


E(o,-L) * Eq(o,-L) + expUKo)] (14) 

where E^(o,-L) is the unperturbed field and 'j'(o) is the random 
phase function. Tl:»e irradiance at a point r due to this source 
field is then from Eqs,(12/j (13) and (14), 


ivr/ 


V V 

is(p,z; J5 KP,z} 


« (k^ I ) [4n^b (z + L)^ (1 + 4L^^)]"^ 
o 


* // dS exp{ik[|p-Si|^ - |p-J2|2]/2(.+l)1 

X expi-ktii^/(b-2iL) - kd2^/0>+2iIi) + iWoi) - iiKoa)! (15) 
and the mean irradiance 


<!(?)> = k^I^[4TT^b(z + L)^(l + 4L^/b^)]“^ 

X // d^Oj d^02 exp{-kOj,^/(b-2iL) ~ k02^/(b + 2iL) 

+ ik[|pA|' “ |p-02p/2(z-^L)l 

X <exp[iiK0l) - i’J^(P 2 )]> 


( 16 ) 
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This differs from th?^. Came quantity in the abaenre of the 
turbulent layer (phase screen) only by the last terra, 

<exp [i'I'C oj) - it|'(02)]>» and its form is determined by the 
characteristics of the turbulent layer* 

The index of refraction in a turbulent region can be 
described in statistical terras, It is usual to write it as 
n(r) •» 1 + nx(r) with ni(r) being a zero mean Gaussian random 
function and to assume that the turbulence is isotropic and 
statistically homogeneous. In this case the two point correlation 
function of the index of refraction fluctuations 


Bjj(ri,r 2 ) = <ni(ri> m(r 2 )> (17) 

depends only on the separation of the two points 


B^Cri-ra) 




ri-r2 ) 


B„(r) 


(IS) 


The spatial spectrum of the index of refraction (fluctuations, the 
3-D Fourier transform of Bj^(r) , is commonly taken to be the 
modified Von Karroan spectrura^^ 

(})„(K) C + K 6xp(-kVk ^). (19) 

n n 0 m 

Here and are the wave numbers corresponding to the 
largest and smallest scale sizes of the turbulent eddies 
respectively and is a measure of the strength of turbulence 
related to the variance <ni > by the fact that 
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( 20 ) 


<nj^> - Bj^(O) - / d^K 

For, if Kg « Kjfl performing the integration yields 

<ni^ - 2Trr(3/2)ra/3)[r(U/6)]"^K^’“^^^ 

- 15.9 K c/ (21) 

o n 

The random ph . screen corresponding to a thin layer of 
turbulence of thickness Az with the spectrum of Bq,(10) is then 
such that^5 

<exp[iK^l) - i'l'(52)]> “ expl-2 [<'|'^> ~ B^|^( oi, 02 ) ]} 

- exp 1-4 A ^Az / Am. *^(K)[l - J^(K 012)]) 

H expl“D^/oi2>} (22) 

1 *^ J 

01-02 > B^(oi, 02 ) is the phase correlation function 

in the plane and 0 ^( 012 ) is called the phase structure function, 

The variance of the phase^^ 

<A “ 2tt^ k^Az / dKK <f> (K) 

n 

» (1271^/5) K G ^ k^ Az 
o n 

« 23.69 K C ^ k^ Az. (23) 

o n 

There are two regimes for which the expression for D((, 
simplifies. If OK^n « 1, then the Bessel function may be 
expanded in series and the integral performed term by term to 
yieldl-B 
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mrnm pm 3 m 
<9f mm mmw 


- 27.48 0^^ 4z - 1.09 °IZ^ (24) 

On the other hand if cKn, » 1 the exponential factor in 

may be replaced by unity and the integral evaluated to give for 

the phase structure function^^*^® 

D^(012) " 44.06 Az 012 ^^^ 

« 1.86 <'|'S 

s (25) 

These two regimes correspond to two quite different physical 
situations, The first corresponds to the case where phase 
correlations are of interest only over a span which is smaller 
than K„j \ i.e., the scale size of the smallest turbulent 
eddies. This would be the case if the diameter of the laser beam 
at the phase screen were smaller than The quadratic 

behavior describes the situation in which the beam is deflected 
but the wavefront across the beam is not distorted and the beam is 
not otherwise affected, The second regime corresponds to the 
case of a beam much larger in diameter than the smallest eddies 

5/3 

and the 012 dependence indicates that the wavefronts across 
the beam are distorted leading to an increase in the diameter of 
the beam farther down the propagation path and a consequent loss 
of power density. 
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The effect that a phase screen or a thin turbulent layer has 
on CARS depends on several factors. For the case of coaxial 
focused Ga’iissian beams with the beam waists located some distance 
from the phase screen as illustrated in Figure 1 both the pump and 
Stokes waves experience identical phas'>. chanties. Hence there is 
no loss of coherence relative to each other, the wave vector 
matching condition is not affected and only the loss in power 
density in the focal region of the beams due to the turbulence 
affects the anti-Stokes power. If the beam diameter at the phase 
screen is smaller than the smallest turbulence scales the only 
effect will be a random steering of the beams without loss of 
power density and the anti-Stokes radiation will be unchanged 
except for a random variation in the propagation direction. On 
the other hand if the laser beams are not coaxial as is the case, 
for example, in BOXCARS, the random deflection of the pump and 
Stokes waves are not identical and may be completely uncorrelated 
depending on their separation and the largest turbulence scale 
size. In this case a considerable disruption of the coherence and 
overlap of the laser beams may occur and the resulting CARS signal 
be severely affected. 

Since deflection of coaxial pump and Stokes beams without a 
change in the instantaneous beam waist diameters does not change 
the magnitude of the measured CARS signal there is no appreciable 
effect in the case where the beam diameters in the turbulent layer 
are smaller than the smallest turbulence scale, Kj„ ^ . In 
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practice this condition would obtain only in very large flow 

tubes. In most cases it is necwo^ary than to consider the effect 

» 

on the pump and Stokes power densities in the focal region 
generated by a phase screen whoso correlations are described by 
Eq.(25). 

Insertion of Eq,(25) into Eq.(16), making the change of 
variables p o « Oj - 02 ^ S » Oj + 02> arid integrating over S gives 

I(r) “ I^k’^[AiT^(z + L)^] ^ / d^o exp {•^iko*r/(z + L) 

- k(4z^ + b^)oV[8b(z + L)^] - 0®/^} (26) 

The angular part of the integration can be performed in turn to 

yield^O 

00 

I(r) = I k^[2TT(z + L)^]"^ / do 0 J [kro/(z + L)] 

o 

X expl~ko^(4z^ + b^)/[sb(z + L)^] - (27) 

Since the CARS signal is produced mainly in the focal region, z » 

2 2 5/3 

0 , and typically kbo /L » A\)jO the turbulence term may 
be regarded as small and exp expanded in series. 

The remaining integration in Eq.(26) may then be accomplished term 
by term. Retaining only the first two terms one has 
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l(r) *• TT ^ I kb(4z^ + b^) ^ {exp [2kbir^/(4z^ + b^) ] 


- A^r(li/6)[k(4z^ + b^)/(8b(z + 

X iFi(ll/6;l; - 2kbr^/(4z^ + b^))} 

® IT ^ I^kb(4z^ + b^) ^ exp[~2kbr^/(4z^ + b^)] 

X (l - A^r(ll/6)[k(4z^ + b^)]”^''® [8b(z + 

X iFi(-5/6;l; 2kbr^/(4z^ + b^))} (28) 

where iFi(a;b;z) is the confluent hypergeoraetric function and 
Rummers transformation^^ has been used. The distribution I(r) 
given in Eq.(28) is nearly gaussian and may be conveniently 
approximated by 


Il(r) = I kb^8“\z) (4z^ + 
o 

X exp{~2kb^r^[3(z)(4z^ + b^) ] (29) 

with 

3(z) = b{l + (5/6) A,r(ll/6) (8b/k)^^® (z + L)^^^(b^ + 4z^)"®'^®) 

(30) 

The factor 3(z) represents additional beam spreading due to the 
phase screen. For example at the beam waist, z = 0, the beam 
radius is larger by a factor [l + (5/6) A^r(ll/6) (8/bk)^'^^ 

The appropriate expression for the field in the presence of a 
phase screen at z = -L instead of that given in Eq.(i) is then 
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E.(?) - (k.I./n)^/2 ^ 2iz/b)“^ 

kJ J iJ 

X exp(ik.z ~ k.p^ b[P(z)(b + 21z)]“^) (31) 

tJ tJ 

and this expression should be used to represent the fields of the 
pump and Stokes beams for calculation of the CARS intensity as 
given in Eq,(7). Although P(z) does depend somewhat on the 
wavelength of the radiation through the k~^^® factor the 
difference between the pump and Stokes wavelengths is sufficiently 
small that one may use the pump wavevector, kj, throughout with 
negligible terror. 

IV . Loss of CARS Intensity Due to a Turbulent Layer 

The total radiated anti-Stokes power in the absence of 
turbulence may be calculated from Eqs.(8), (9), and (10). The 
effect of a turbulent layer in the path of the pump and Stokes 
driving waves remote from the focal region is to induce a 
reduction in the on-axis amplitude and an increase in the 
effective diameter of the driving beams. Substitution of the 
expression given in Eq.(31) for the pump and Stokes fields in 
Eq.(7) and evaluation of the integrals over < and p yields an 
equation similar to Eq.(8); 

E3(r^) = i(3NXk^^)Tr"’^/2 (4k3)“l(kjj exp(ik3Z^) 

z 

o , 

X / dz(l + 2iz/b)“^ $“^'^(z)[2ki +k 2 - 2 ik 3 z/b]"^ 

-L 

^ exp[-p^^/bH^(z^,z) ] (32) 
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with 




4z^/b^) 


b(2ki 


g(z )(l -f 

+ kz “ 2ik3z/b) 


- 2i 


(z “ 




bk3 


(33) 


The total anti-Stokes power radiated is then 

l3(Lg) « 27T 7 dPoPo|E3(P^;V L^)|^ (34) 


which may be rewritten using Eq.(32) as a double integral 
l3(Lg) ** gN^X^k^*^ (4Tik3)’‘^ ki^I 1^212 


X / / dzidzad + 2izi/b) Hi - lizz/h) ^ 

-L 

X S”^''Hzi) 3 “^'^^(z 2) [2ki + k2 - 2ik3Zi/b]"^ 

X [2ki + k2 + 2ik2/b]“- [h^(L^,zi> + H^*(L^, Z 2 )]”" (35) 


As stated earlier at the end of Section II the major contribution 
to the CARS signal arises within a distance equal to ten times the 
confocal parameter of the beam waist. All lengths along the 
propagation path may be scaled in units of the confocal parameter, 
i.e,, the dimensionless variables xi,2“ 2 zi, 2 /b may be defined 
and the range of integrations restricted to -10 < xi,2 < 10 to 
give assuming L > 5b and Lg > 5b, 


ORI01NAL pmn m 
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I3CA ) * 9 NVk ki^ k2 I2 b 

e o 


« dxidx2 (1 + ixi)”^(l - iX2)"V^^^ (xi) y’^^^(x2) 
-10 


X [ 2 ki + k 2 - ik 3 Xi]~^[ 2 ki + k 2 + ik3X2] ^ 

X [h^(A^,xi) + (36) 


where H ■ 2L /b, A » 2L/b, 

e e , 

•r(x) “ 3(x)/b «= [l+(5/6)A^r(ll/6)(2b/ki)"'^/®(x+^)®/^(l+x)'"^''®] 


and H^(A^,x) 


Y(x)(l+x^) _ . ^ 

2ki+k2-2ik3X ka 


(37) 


The integrals in Eq.(36) must be evaluated numerically. 
However, if H » 10, i.e. L » 5b, within the range of integration 

Y(x) « Y(o) = 1 + (5/6)A^r(ll/6)(2b/k)^^® ^ (38) 

“I /2 “* 1/2 

Moreover the Y (xi)Y (x 2 ) factors are dominant compared to 
those involved in Hj, and . In this approximation the signal 
in the presence of a turbulent layer is just Y ^(o) times that in 
the absence of turbulence, i.e., the ratio 
I 3 /I 3 = Y(o) 

= 1 + (5/6)A^r(ll/6)(2b/ki)^^® 

= 1 + (5/6)A^r(U/6)(8/bki)^''® 



Thus, in this approximation, the ratio of the CARS signal in the 
absence of turbulence to that generated through a turbulent layer 
varies linearly with the strength of turbulence through the factor 
A \|>5 it varies as b ' so that tighter focusing, i.e., 
decreasing b, is detrimental; and the 5/3 power dependence on the 
distance from the turbulent layer to the focal region imposes a 
restraint on the region which may be explored through a turbulent 
layer. The validity of the approximation and the behavior 
expected for realistic estimates of the parameters 
involved is illustrated in Figures 2 and 3 where ^> 0 ( 13 / 13 ) 
calculated from Eq.(39) and the value determined from numerical 
integration of Eq.(36) are plotted as functions of the distance 
between the phase screen and focal region, L. The pump, Stokes 
and anti-Stokes wavelengths are assumed to be 532 nm, 557 nm and 
509 nm respectively for purposes of this illustration. 

The values of A\jj used to generate the data displayed in the 

figures were calculated on the basis of some assumptions of the 

characteristics of the turbulent boundary layer. Referring to 

Eq.(25), A)|; = 44.06 ki^ Az, Az being the thickness of the 

2 

turbulent layer and being the strength of turbulence 
parameter which is related to the variance in index of refraction 
<ni > and the wave number of the largest scale turbulent eddies 
Kf, through Eq.(21). For this illustration it is assumed that 
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the phase screen models the turbulent boundary layer in a mach 2,5 
flow in a rectangular nozzle of 50 cm cross-sectional area. 

k 

Typical characteristics for this type of flow^S are 

Az * 5 mm, K '*3,1 cm”^, K * 630 cm" ^ and the rms density 
o m 

2 1/2 

fluctuations <(Ap/p) > may be as great as 5%, The index of 

refraction of air varies linearly with density with the deviation 

from unity, ni “ 3 x 10 ** (Ap/p). Thus <niS < 2 x 10 

2 -11 -2/3 

the strength of turbulence parameter < 2 x ],0 cm , the 

2 « c/ 3 

phase variance <i|^ > < 0.5 and < 6 cm . The values of used 

— 5/3 . 

in the following calculations are 0.6 and 6 cm corresponding 
to rms density fluctuations of 0,5% and 5%. 

Two values of the confocal parameter were also chosen for 


purposes of lllu'^trat lont 


mm and 0.14 mm, corresponding to 


focusing a 1 cm diameter collimated beam with 40 and 10 cm focal 
length lenses respectively. Rewriting Eq,(39) as 
I3/I3 B 1 + (40) 

with o = 4,43 A^(bki) the value b = 0.14 mm generates values 
of a = 5.5 X 10 "^ and 5.5 x io~^ cra“®^^ for A^ “ 0.6 and 6 
respectively. Curves of £ 0 ( 13 / 13 ) calculated from Eq.(40) with 
these values of o are plotted in Figure 2 along with the numerical 
integration results. Similar curves for b = 2.2 mm, 
a = 5.5 X 10 ** and 5.5 x 10 ^ cm ^ are plotted in Figure 3. In 
both cases it is seen that a significant loss of signal can 
occur. For example, 'rms density fluctuations of 5% cause a factor 
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of 2 loss in signal if the turbulent layer is 5.7 cm from the 
focus of a 10 cm focal length lens (Figure 2). 

V. Discussion 

It is clear from Figures 2 and 3 that the approximation 
invoked in deducing Eq,(39) from Eq,(36) is reasonably accurate 
over a substantial range of separation of the phase screen and 
focal region, L, turbulence strength, Ajj,, and beam geometry, b. 
Based on the indicated functional dependence it may be coniiluded 
that it is advantageous to reduce as far as possible the 
separation of the focal region and the turbulent layer and to 
focus the beam less tightly, i.e,, increase b. Of course it is in 
gtneral not possible to choose the distance between the focus and 
the boundary layer. Nor is it always desirable to increase the 
confocal parameter since that implies an increase in the beam 
diameters and the length of the region in which the signal is 
generated with a concomitant loss in spatial resolu<’ion. 

The results presented were calculated for turbulence 

strengths and scale sizes which are typical of transonic flows in 

small nozzles and the tightness of focus of the driving beams was 

chosen to conform to a typical laboratory experiment. Application 

of these results to CARS in flames, jets and other high speed 

flows should also be possible. In order to do so one would need 

2 , l/ 2 

values for the rms density fluctuations, <(/\p/p) > ' , in the 
turbulent boundary between the flame or jet and the surrounding 
atmosphere as well as the thickness of the layer, Az, and the 
largest turbulence scale size, . The index of refraction 
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A 

variance <ni can be calculated if dn/dp is known and then 
Eqf(21) and (25) used to find 

*• 1,5 <ni^> ^ Az from which « ■ 4.4 A^(bk) 

can be determined. The ratio of the anti-Stokes power in the 

presence of turbulence to that without is then 
*• 5/3 —1 

1 3/1 3 ■ (1 + otL ' ) with L being the distance frotn the beam 
waists to the turbulent layer. This simple calculation should 
indicate whether or not a CARS experiment is feasible in a given 
situation. 

It is also v/orth reiterating that if the beam diameters as 
they penetrate the turbulent layer are smaller than the smallest 
scale phase perturbations then the beams are simply deflected 
without loss of instantaneous power density or coherence. In the 
coaxial beam case the CARS process is completely unaffected except 
for a small fluctuation in the direction of the emerging 
anti-Stokes radiation. If| however, a crossed beam geometry is 
used the pump and Stokes beams are deflected independently with 
the result that the average overlap of the beams is reduced and 
the signal degraded. 
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Figure 1. Illustration of the increase in the mean diameter of a 
beam wave due to a turbulent layer. 
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Figure 2, Logarithm of the ratio of the unperturbed anti-Stokes 

power to that perturbed by a phase screen versus 
5/3 

L . A, 0 - values calculated by numerical 
integration of Eq.(36). Solid lines: An(l + 

with a = 5.5 x lO'"^ and 5.5 x lO"^ cra"^^^ 
corresponding to rras density fluctuations of 0.5 and 
5% respectively. Confocal parameter b = 0.14 nin, 
i.e., 1 cm diameter beam focused by a 10 cm fl lens. 
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Figure 3, Same as figure 2 except b •= 2.2 ram corresponding to a 
1 cm diameter beam focused with a 40 cm fl lens, In 
this case o. ® 5.5 x 10 ** and 5.5 x 10 ^ for 
<(Ap/p)^>^^^ « 0.5 and 5% respectively. 
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